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1. Introduction 

When negative muons pass through matter, they can be captured into high-lying 
atomic orbitals. Then, in a time-scale of the order of 10~ 13 s, they cascade down 
into the Is orbit, through Auger processes with atomic electrons and the emission 
of X-rays. At this point, two competing processes occur: one is ordinary decay 

(j," ->• e~ + V e + , (1) 

and the other is the (weak) capture by the nucleus 

p- + A{Z,N)^v„ + A(Z-1,N + 1) , (2) 

which can take place from any of the two initial hyperfine states, / = Jj ±1/2 (Jj is 
the spin of the initial nucleus A(Z,N)). Apart from tiny corrections due to bound- 
state effects ffl, the decay rate is essentially the same as for a free muon. In light 
nuclei, this is much larger than the rate for capture, which proceeds predominantly 
through the basic process 

M~ + p -> n + , (3) 

induced by the exchange of a W + boson. Its rate is expected to be proportional to 
the number of protons in the nucleus and to the probability of finding the muon 
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at the nucleus. Since the semi-leptonic weak nuclear interaction is effectively a 
contact interaction, this probability scales like the square of the atomic Is wave 
function evaluated at the origin^, proportional to Z 3 . The capture rate, therefore, 
scales roughly like Z 4 . It is only for nuclei with Z > 12 that the nuclear capture 
rate becomes comparable with the decay rate. Muonic capture on light nuclei are 
therefore experimentally challenging processes. However, they are preferred under 
the theoretical point of view, as the nuclear effects can be easier and more accurately 
taken into account, and informations on the basic process of Eq. ([3]) can be better 
extracted. For instance, muon capture on hydrogen and hydrogen isotopes is, in 
principle, best suited to obtain informations on the matrix element of the (charge- 
changing) single-nucleon weak current 

f = «p" 

but is experimentally the hardest process. Note that in Eq. (|4]) we have ignored 
contributions from second-class currents for which there is presently no firm 
experimental evidence^, and we have indicated with u p (u n ) the proton (neutron) 
spinor, with Mjy the nucleon mass and with q 2 the four-momentum transfer; 7 M 
and 7 5 arc the standard Dirac matrices, a^ iV = 1/2^, j v ] and — u^jo • Of 
the four form factors of Eq. (j4}, Fi(q 2 ) and i*2(<Z 2 ) are related to the isovector 
electromagnetic form factors of the nucleon by the conserved- vector-current (CVC) 
constraint. They are well known over a wide range of momentum transfers q 2 from 
elastic electron scattering data on the nucleon^. The axial form factor Ga(q 2 ) is 
also quite well known: its value at vanishing q 2 , qa = 1.2695 ± 0.0029, is from 
neutron /3-decay^, while its (^-dependence is parametrized as 

G A (q 2 )=g A /(l-q 2 /A 2 A ) 2 , (5) 

with A A = 1 GeV from an analysis of pion electro-production data ^ and direct 
measurements of the reaction p + — > n + /i + ®. Note that a considerably larger 
value A A — 1-35 GeV is obtained from current analyses of neutrino quasi-elastic 
scattering data on nuclear targets However, these analyses are based on rather 
crude models of nuclear structure (Fermi gas or a local density approximation of 
the nuclear matter spectral function) and on simplistic treatments of the reaction 
mechanism. Also, some discrepancies exist on the neutron /3-decay lifetime, as the 
world average value used here ^differs by 6.5 standard deviations from the results 
obtained from gravitationally trapped ultra-cold neutrons H[ A discussion of this 
point is however well beyond the subject of the present review. 

The induced pseudoscalar form factor Gps(q 2 ) is the least known of the four 
form factors of Eq. ((4]). The MuCap collaboration at Paul Scherrer Institute (PSI) 
has recently reported a precise measurement of the rate for reaction Q in the 
singlet hyperfme state (/ = 0): 725.0 ± 13.7(stat) ± 10.7(syst) sec -1 ^A. Based 
on this value, an indirect "experimental" determination of Gps at the momen tum 
transfer q$ = — 0.88m^ relevant for muon capture on hydrogen has been given 
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Gpg P ((/Q) = 7.3 ± 1.2, by using for the remaining form factors the values discussed 
above and by evaluating electroweak radiative corrections. These are found to be 
sizable, of the order of ~ 3 %. Theoretical predictions for the induced pseudoscalar 
form factor were derived long ago based on the notion of a partially conserved 
axial current (PCAC) and pion dominance, and were later refined by evalua ting 
leading-order corrections to the PCAC result with current algebra techniques ^1 
More r ecently , these predictions have been re-derived in chiral perturbation theory 
(XPT) [ISEST finding 

GpsiQ ) = 2 2 ^9Am a M N r A , (6) 

where g wpn = 13.05 ± 0.20 is the ttNN coupling constant, /„. = 92.4 ± 0.4 MeV is 
the pion decay constant, and r A — 0.43 ±0.03 fm 2 is the axial radius of the nucleus, 
related to A A of Eq. © as = 12/r|. For ql = -0.88 m 2 , G™(g 2 ) = 8.2±0.2E3 
To be noticed that the evaluation of electroweak radiative corrections for muon 
capture on hydrogen is crucial for bringing Gp* p within less than la of G™. 

Besides their relevance for extracting informations on single-nucleon weak cur- 
rent form factors, muon captures on light nuclei also provide a testing ground for the 
theoretical frameworks used to study those reactions of astrophysical interest whose 
rates cann ot b e measured experimentally, and for which one has to rely exclusively 
on theory 03. In fact, the same nuclear wave functions and, indirectly, the same 
model for the nuclear interactions from which these are obtained, and the sam e nu- 
clear weak current can be used to study neutrino reactions in light nuclei weak 
proton captures on proton and 3 He (the so-called pp and hep reactions) ' ' ' \ 
and muon captures on light nuclei. 

In the present review, we focus our attention on the following captures: 

fx" + d -> n + n + , (7) 
pT + 3 Hc -> 3 H + i/ M . (8) 

Muon capture on 3 He can also occur through the two- (nd) and three-body (nnp) 
breakup channels of 3 H. However, the branching ratios of these two processes are 20 
% and 10 %, respectively, and experimental and theoretical work on them is quite 
limited. They will not be discussed here. A comprehensive and detailed description 
of these reactions, and more in general of the physics of muon capture and the 
problem of the induced pseudoscalar form factor, can be found in Refs. [22j l23l 
and El 

The observables of interest for muon capture reactions are the capture rates. In 
reaction ([7]) , the stopped muons can in principle be captured from the two hyperfine 
states, / = 1/2 or 3/2. However, it is kno wn tha t capture takes place practically 
uniquely from the doublet hyperfine state I 22 | 23 | therefore consider only 

the doublet capture rate T D . In reaction (JHJ, instead, a difference in the capture 
rates between the hyperfine states is not expected: a hyperfine transition is highly 
unlikely, due to the energy difference between the hyperfine states. We will therefore 
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consider the total capture rate Tq. In the next section we will briefly discuss the 
experimental situation, while the formalism to derive these observables and the 
most recent theoretical calculations will be presented in Sec. [3] Some concluding 
remarks are given in Sec. [4] 

2. Experimental situation 

T he fi rst attempt to measure T D was carried out over forty years ago by Wang et 
aL^Sl. Using a liquid mixed H2/D2 target, they obtained T D = 365 ±96 s" 1 . A few 
years later, Bertin et al. measured T D = 445 ± 60 s _1 t2S] using a gas mixed H 2 /D 2 
target, and assuming a pure doublet mix of [id spin states. However, a subsequent 
study of hyperfine depopulation in a H 2 /D 2 mixture ^ has failed to support this 
assumption, and therefore the Bertin et al. result is considered controversial^!. The 
most recent measur em ents have been performed in the eighties by Bardin et al. ^Sl, 
and Cargnelli et al. ^1 They both used pure deuterium, so that the /id + d collision 
rate is sufficient to fully depopulate the / = 3/2 hyperfine state. Furthermore, 
Bardin et al. used a liquid target and the lifetime method, i.e. compared the negative 
and positive muon lifetime when stopped in deuterium. For positive muons, the 
lifetime is the inverse of the muon decay rate, while for negative muons the lifetime 
is the inverse of the sum of the muon decay rate and the muon capture rate. The 
lifetime difference thus determines the capture rate, assuming, according to the 
CPT theorem, the positive and negative muon decay rate identical. Cargnelli et 
at, instead, used a gas target and the neutron method, i.e. directly detected the 
recoil neutrons, which is obviously quite challenging. These two measurements gave 
T D = 470±29 s _1 l2Hland r D = 409±40 s -1 ^. In conclusion, all the measurements 
available until now, while consistent with each other, are not very precise, since 
the errors are in the 6-^10 % range. Howev er, there is hope to have this situation 
clarified by the MuSun Collaboration! 3 -^, with their on-going experiment at PSI, 
which should reach a precision of 1.5 %. The gain in experimental precision relies 
on the fundamental techniques developed for the MuCap experiment Muons 
will be stopped in an active gas target consisting of a cryogenic ionization chamber 
operated as time projection chamber with ultra-pure deuterium. The muon stopping 
point will be reconstructed in 3 dimensions, and this will eliminate the otherwise 
overwhelming background from muon stops in wall materials. The capture rate will 
then be determined using the lifetime technique. 

The experimental situation for muon capture on 3 He is much clearer. After a 
first set of measurements in the early sixties by Falomkin et al. ^H, Zaimidoroga 
et al. Auerbach et al. ^221 anc | Clay et al. a very precise determination was 

Iq r I 

performed by Ackerbauer et al. in the late nineties. The basic method involves 
counting the numbers of muon stops and 3 He recoils, when a beam of muons is 
stopped in 3 He. Ackerbauer et al. used a gas ionization chamber, which allowed a 
better separation of muon and 3 He signals. The measured total capture rate To, 
corresponding to a statistical average population of the 4 different hyperfine states, 
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is 1496 ± 4 s _1 1"*", a value consistent with those of the earlier measurements, but 
with a factor 10 of improvement in the experimental accuracy. 

If the hypcrfine structure of the (/i, 3 He) system is taken into account and the 
direction of the recoiling triton is detected, there are, in addition to the total cap- 
ture rate, other observables, i.e. angular correlation parameters or so-called recoil 
asymmetries, which are more sensitive than the capture rate itself to the value of 
the induced pseudoscalar form factor Gps(q 2 )- A first attem pt t o measure the re- 
coil asymmetry has been made by Souder et al. at TRIUMF l^fil. They used a 3 He 
ionization chamber to stop the incoming muons, re-polarize the (/i, 3 He) system and 
track the triton recoils. They obtained for the vector asymmetry A v the value of 
0.63 ± 0.09 (stat.)^Q 14 (syst.). This experimental result, which to our knowledge 
represents the first measurement of this observable, is affected by large systematic 
uncertainties. Therefore, a comparison between theory and experiment would not be 
particularly meaningful. Thus, further experimental work is highly recommended. 



3. Theoretical calculations 

Before discussing the results of the different theoretical calculations for the capture 
rates of reactions (J7|) and (O, we present in the following subsection the formalism 
used in the calculation of the observables under consideration. 



3.1. Theoretical formalism 

The m uon capture on deuteron and 3 He is induced by the weak interaction Hamil- 
tonian S3 

H w = < ^Jtel <J {*)f{*), (9) 

where Gy is the Fermi coupling constant, Gy=l. 14939 xl0~ 5 GeV~ 2 as obtained 
from an analysis of + — > + /3-decays an d l„ and j a are the leptonic and 
hadronic current densities, respectively. The former is given by 

Z CT (x) = e - ik "- x u(k l/ ,/i l ,)7 CT (l -7 5 )Vv( x : s a0 > ( 10 ) 

where Vv( x > s m) * s ^ ne ground-state wave function of the muon in the Coulomb field 
of the nucleus in the initial state, and uQz„, h v ) is the spinor of a muon neutrino with 
momentum k y , energy E v (=fc„), and helicity h v . While in principle the relativistic 
solution of the Dirac equation could be used, in practice it suffices to approximate 

k„-X) , (11) 

since the muon velocity ~ Za <C 1 (a is the fine-structure constant and Z=l or 2 
for deuteron or 3 He, respectively). Here ipi s (x) is the Is solution of the Schrodinger 
equation and, since the muon is essentially at rest, it is justified to replace the two- 
component spin state x( s m) w hh the four-component spinor u(k^,s^) in the limit 
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k,, — > 0. This will allow us to use standard techniques to carry out the spin sum 
over s p at a later stage. 

In order to account for the hyperfine structure in the initial system, the muon 
and deuteron or 3 He spins are coupled to states with total spin /, equal to 1/2 or 
3/2 in the deuteron case, and to or 1 in the 3 He case. The transition amplitude 
can then be conveniently written as ^21 

Tw(f,fz;si,s 2 ,h v ) = (nn,s 1 ,s 2 ;v,h l/ \H w | (p,d);f,f z ) 

(* PlS1S2 (nn)|r(q)|* d ( Sd )) , (12) 
for the muon capture on deuteron, where p is the nn relative momentum, and SOI 
T w (f, f z ; 4, K) = ( 3 H, 4; u, K\H W \ ( M , 3 He); /, /,) 

(*3 H (4)li CT (q)I^He(s 3 )) , (13) 

for muon capture on 3 He. In Eqs. (|12[) and (fT3| we have defined 

la(K, s^) = u(k„, h v ) j a (1 - 7 5 )u(k M , s M ) , (14) 
and the Fourier transform of the nuclear weak current has been introduced as 

j CT (q) = J dxe^ x r(x) = (p(q)J(q)) , (15) 

with the leptonic momentum transfer q defined as q = k^ — k„ ~ — k^. The function 
ipis (x) has been factored out from the matrix element of j CT (q) between the initial 
and final states. For muon capture on deuteron, ipfl is approximated as ^3 

IKI 2 - hM0)l 2 = (16) 

where V'is(O) denotes the Bohr wave function for a point charge e evaluated at the 
origin, and fi^d is the reduced mass of the (//, d) system. For muon capture on 3 He, 
ipf s is approximated as^^l 



ICI 2 ^ ^ , (17) 

7T 

where in this case n^ne is the reduced mass of the (/i, 3 He) system, and the factor 
TZ approximately accounts for the finite extent of the nuclear charge distribution I27J. 
This factor is defined as 



^ = U. 7nM2 ' ( 18 ) 

with 
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where p(x) is the 3 He charge density. It has been calculated explicitly in Ref . l39l bv 
using the charge densities corresponding to two realistic Hamiltonian models, the 
AV18/UIX and N3LO/N2LO (see below), and has been found for both models to be 
within a percent of 0.98, the value obta ined from the experimental charge density 
and commonly adopted in the literature ^3. 

In the case of muon capture on deuteron, the final state wave function is ex- 
panded in partial waves as 



* P ,. 1 ,. a (nn)=4w^(±a 1 ,±a 2 \SS,) ]T i L Y^ (p){SS z , LL Z \ JJ Z ) <f J ' (p) , 

S LJJ, 

(20) 

XjSJJ 

where ^> nn z (p) is the nn wave function. The calculation is typically restricted to 
J < 2 and L < 3, sin ce it has been proven that higher order partial waves give 
negligible contributions^. Therefore, in spectroscopic notation, only the 1 Sq, 3 Pq, 
3 Pi, 3 P2- 3 F 2 , and 1 D 2 parti al wave s are considered. 

Now, standard techniques | 20 | 37 | are use d to carry out the multipole expansion of 
the weak charge, p(q), and current, j(q), operators. For muon capture on deuteron, 
we find 



(Kn JJ *(p)\p(cL)m(s d )) = ^Y. i A ilSd ^ z) c^( q ) , (2i) 

A>0 V2J + 1 

<^2 JJ- (p)|j.(q)l*d(-)> = -V4^E ^A^l i A (1Sd ^ z) Li SJ ( q ) , (22) 

A>0 V2J + 1 

A>1 V <W + 1 

[-XM^ J (q) + E^ J (q)} , (23) 



where A = ±1, and C^ SJ (q), L^ SJ (q), E^ SJ (q) and Mj( SJ (q) denote the reduced 
matrix elements (RME's) of the Coulomb (C), longitudinal (L), transverse elec- 
tric (E) and transverse magnetic (M) multipole operators, as defined in Ref. [20] 
Since the weak charge/current operators have scalar /polar- vector (V) and pseudo- 
scalar/axial- vector (^4) components, each multipole consists of the sum of V and 
A terms, having opposite parity under space inversion 1201. The contributing multi- 
poles for the S-, P-, and £)-channels mentioned above in muon capture on deuteron 
are given in Table [1] where the superscripts LSJ have been dropped. In the case 
of mu on capture on 3 He, explicit expressions for the multipole operators are given 
bySQl 



(*8 H (4)|p(q)l*»Ho(*3)> = V^E ^TTii 1 <4,o(-0) (^3,/m|^4)Q(?) 

z=o,i 



(24) 



December 2, 2011 1:34 WSPC/INSTRUCTION FILE mucap-review 



8 Laura Elisa Marcucci 



(*»H(4)li.(q)l*»He(* 3 )) = -V&F E VWTTi 1 d l mfi (-e) (ls 3 , i rn\l s ' 3 ) lm) , 



(*. H (*3)|jA(q)|*»He(*3)) = VS^ i<_ A (-0) (-*3,lm|-4)[-AMi( ? ) + £?i(g)] , 

(26) 

where m=s 3 — S3, and the d l mm , are rotation matrices in the standard notation 
of Re f. 1411 Applying parity and angular momentum selection rules, it has been 
shown 40 that the only contributing RME's are C (V), Ci(A), L (V), L X (A), E^A), 
and Mi(V). 

The total capture rate for the two reactions under consideration is then defined 

as 

dT = 2n5(AE)\T w \ 2 x (phase space) , (27) 

where 5(AE) is the energy-conserving (5-function, and the phase space is 
dp dk„/(27r) 6 for reaction (J7J and just dk t/ /(27r) 3 for reaction (J8)). The following 
notation has been introduced: (i) for muon capture on deuteron 

w = -ttt E £i2v(/>/*;«i,«2,M a , (28) 

S1S2/I1/ /a 

and the initial hyperfine state has been fixed to be / = 1/2; (ii) for muon capture 
on 3 He 

W = EE p (/' /»; 4. Ml 2 , (29) 

s^/l„ //a 

where P(/, / z ) is the probability of finding the (/x, 3 He) system in the total-spin 
state f,f z and P(f, f z ) = 1/4 when the same probability to the different hyperfine 
states is assigned. 



Table 1. Contributing multipoles in muon cap- 
ture on deuteron, for all the nn partial waves 
with J < 2 and L < 3. The spectroscopic nota- 
tion is used. See text for further explanations. 

Partial wave Contributing multipoles 



1 S Ci(A), Li (A), E-i(A), M 1 {V) 

3 P Ci(V), Li(V), Ei 00, Mi (A) 

3 Pi C (A), L (A), 

Oi(V), Li(V), Ei(V), Mi(A), 
C 2 (A), L a (A), E 2 (A), M 2 (V) 
3 P 2 - 3 F 2 Gi(X0, LiOO- E i(V), Mx (A), 
C 2 (A), ^2(A), £ 2 (A), M 2 (y), 
C 3 (V), L 3 (V), E 3 (V), M 3 (A) 
1 D 2 Gi(A), Li(A), Ei(A), M 1 {V), 

C 2 (V), L 2 {V), E 2 (V), M 2 (A), 
C 3 (A), L 3 (A), S3 (A), M 3 (V) 
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After carrying out the spin sums, the total rate and recoil asymmetry for muon 



40 



capture on 3 He are 

T = Gy E v ( 1 - 



mi 



\C (V) - L (V)\ 2 + \d(A) - L^A)] 2 + \Mi(V) - E 1 (A)\ 2 , (30) 
with E v given by 

(rrift + rn3 Hc ) 2 - m 2 H 



E„ 



2(m p + m3 Ho ) 



and 



21m 



.4 



1 



(C (V) - LoOTftfaiA) - - \Mx{V) - E 1 (A)\ 



\C (V) - L (V)\2 + \Ci(A) - LMW + \Mi(V) - E 1 (A)\^ 
In the case of muon capture on deuteron, the differential rate reads 



(II 



dp 



° =El 



1 - 



E v 



mi 



,p 2 dp 



where 



E„ 



(m M + m d ) 
M + m d ) 2 - Am 2 n - 4p 



8tt 4 



(31) 



(32) 



(33) 



(34) 



2(m M + m d ) 

In Eqs. (|30 [i ~ (|34|) . m^, m n , m d , ?Ti3Hj ™shc are the muon, neutron, deuteron, 3 H 
and 3 He masses. The integration over p in Eq. ((33 j) is performed numerically using 
Gauss-Legendre points. A limited number of them, of the order of 10, is necessary 
to achieve convergence to better than 1 part in 10 3 . In order to calculate the total 
capture rate T D , the differential capture rate is plotted versus p, and numerically 
integrated. Usu ally, about 30 points in p are enough for this integration in each 
partial wave ^21. 



3.2. Results 

Theoretical work on reactions ([7]) and © is just as extensive as the experimental 
one (see Sec. [2]). A list of publications, updated to the late nineties, is given in 
Table 4.1 of Ref. [22j in Ref. [23] and Ref. [24] Here, we limit our considerations to the 
calculations performed since the year 2000. The starting point will be our studies 
of Refs. l39l and 142] for both muon capture reactions under consideration, as, to our 
knowledge, are the most recent on es published. T hese results will be compared w ith 
the calculations of Ando et al. S3 anc j Rj cc i e £ a /. HH for reaction ([7]) , and Gazit ^Sl 
for reaction (jSJ . We will comment also on our early study of reaction ([8]) ^4 and 
on the results of Ho et a/PSl and Chen et aiP^. The most recent studies (available 
only as preprint) of Refs. l48l and l49l will be also briefly discussed. 

The theoretical results of Refs. EU [43] and [44] for the capture rate T D of reac- 
tion (0 from the initial doublet hyperhne state and of Refs. [39] [40] and [45] for the 
total capture rate Tq of reaction ([8]) are summarized in Tables [2] and [3] respectively. 
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Table 2. Summary of the theoretical results 
for the doublet capture rate V D (in s _1 ) of 
muon capture on deuteron. Only the calcula- 
tions after the year 2000 are considered. The 
results obtained with the 1 So nn final state are 
also shown. There where possible, the theoret- 
ical uncertainty is also indicated. 



Rcf. 


r D 


r D Cs ) 


Ando et al. 331 
Ricci et al. S3 
Marcucci et al. EH 


386 
423 ± 7 
392 ± 2.3 


254 ± 1 
261 ± 7 
248.6 ± 2.7 



Let us review our work of Ref. 1391 The first ingredient for any theoretical study 
of the reactions under consideration is the realistic Hamiltonian model used to de- 
scribe the initial and final A = 2 and 3 nuclear wave functions entering in Eqs. (|12[) 
and (I13p . Two representative two- nucleon interaction models have been used, the 
phenomenological Argonne v±s (AV18) EOI anc | the potential derived within chiral 
effective field theory (xEFT) up to next-to-next-to-next-to leading order (N3LO) 
by Entem and Machleidt^l These two models both reproduce the deuteron observ- 
ables and the large two- nucleon scattering database with a \ 2 /datum ~ 1. Given 
the significant differences in their derivation and structure, they are believed to be a 
representative subset of the accurate two- nucleon interaction models available in the 
literature. To accurately describe the A — 3 nuclear systems, it is well known that 
the two-nucleon p oten tials need to be augmented by three-nucleon interactions. The 
Urbana IX (UIX) ^2] model has been used in conjunction with the AV18, and the 
chiral three-nucleon interaction, derived up to next-to-next-to leading order (N2LO) 
in Ref. [53j has been used together with the N3LO. The hyperspherical-harmonics 
(HH) method has been used to solve the A-body bound and scattering problem, also 
in the context of A = 2 systems, for which of course wave functions could have been 
obtained by direct solution of the Schrodinger equation. The HH method for A > 3 
has been reviewed in considerable detail in a series of recent publications I 54 | 55 | 56 l 

The weak current consists of polar- and axial-vector components, derived within 
two different frameworks, the "Standard Nuclear Physics Approach" (SNPA) and 
xEFT. The first one goes beyond the impulse approximation, by including meson- 



Table 3. Same as Table [2] but 
for the total capture rate To of 
muon capture on 3 He. 



Ref. 




r 


Marcucci et al. 


\M 


1484 ± 8 


GazitESl 


ED 


1499 ± 16 


Marcucci et al. 


1484 ± 13 
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exchange currents (MEC's) and terms arising from the excitation of A- isobar de- 
grees of freedom. The second approach includes two-body contributions derived in 
heavy-baryo n chira l perturbation theory (HBxPT) within a systematic expansion, 
up to N3LO ES53 Since the transition operator matrix elements are calculated 
using phenomenological wave functions, this second approach is a "hybrid" %EFT 
approach (xEFT*). Here we briefly review the main characteristics of the weak cur- 
rent operator, both within SNPA and xEFT*. We consider only the contributions 
beyond the one-body term, as the one-body operators can be easily obtained per- 
forming a non-relativistic reduction of the single - nucleon weak current of Eq. (j4|) , 
retaining corrections up to order (q 2 /M^) I 39 | 20 |_ 

The polar (scalar) weak current (charge) operator is related to the isovector part 
of the electromagnetic current (charge) via the CVC hypothesis. In SNPA, no free 
parameters are present in the model for the electromagnetic operator , w hich is able 
to reproduce the trinucleon magnetic moment s to bett er than 1 % "21, as well as 
a large variety of electromagnetic observables I 58 | 59 160| j n cage f ^-eft*, no 
two-body contributions to the scalar charge operator are present at N3LO, while 
the vector current is decomposed into four terms ^H. the soft one-pion exchange 
(l7r) term, vertex corrections to the one-pion exchange (lvrC), the two-pion ex- 
change (2tt), and a contact-term contribution. Their explicit expressions can be 
found in Ref. [57] All the lir, lirC and 2tt contributions contain low-energy constants 
(LEC's) estimated using resonance saturation arguments, and Yukawa functions ob- 
tained by performing the Fourier transform from momentum- to coordinate-space 
with a Gaussian regulator characterized by a cutoff A. This cutoff determines the 
momentum scale belo w which these ^EFT currents are expected to be valid, i.e. 
A=500-^800 MeV The contact-term electromagnetic contribution is given as 
sum of two terms, isoscalar and isovector, each one with a LEC in front (g^s and 
gw), fixed to reproduce the experimental values of A = 3 magnetic moments. The 
resulting LEC's are given in Table V of Rcf. W\ and listed again in Table [H for 
completeness. The uncertainties on (745 and gw are not due to the experimental 
errors on the triton and 3 He magnetic moments, which are in fact negligible, rather 
to numerics. 

Table 4. The LEC's 345 and g^y associated with the isoscalar and 
isovector contact terms in the electromagnetic current, and the LEC 
dn of the two-body axial-current contact term, calculated for three 
values of the cutoff A with triton and 3 He wave functions obtained 
from the AV18/UIX model. For A = 600 MeV, the N3LO/N2LO 
model is also used. 





A (MeV) 


94S 


94V 


da 




500 


0.69±0.01 


2.065±0.006 


0.97±0.07 


AV18/UIX 


600 


0.55±0.01 


0.793±0.006 


1.75±0.08 




800 


0.25±0.02 


-1.07±0.01 


3.89±0.10 


N3LO/N2LO 


600 


O.lliO.01 


3.124±0.006 


1.00±0.09 
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The two-body axial current operators inSNPA as used in Ref. [39l as well as 
in the studies of the pp and hep reactions | 19 | 20 |^ can ^ e divided m t w0 classes: the 
operators of the first class are derived from it- and p-meson exchanges and the pir- 
transition mechanism. These mesonic operators give rather small contributions 
The operators in the second class are those that g ive the largest two-body contri- 
butions, and are due to A-isobar excitation I I I In particular, in the dominant 
A-to-A-transition axial current, the A-to-A axial coupling constant (g A ) is retained 
as a parameter and is determined by fitting the experimental Gamow- Teller matrix 
element of tritium /3-decay (GT EXP ). Also the pseudoscalar term in the A-to-A- 
transition axial current is retained. It is important to note that the value of g* A 
depends on how the A-isobar degrees of freedom are treated. In the muon cap- 
ture studies presented here, the two-body A-excitation axial operator is derived in 
the static A approximation, using first-order perturbation theory. This approach is 
considerably simpler than that adopted in Ref. [20l where the A degrees of freedom 
were treated non-perturbatively, within the so-called transition-correlatio n o pera- 
tor approach, by retaining them explicitly in the nuclear wave functions ^B. The 
results for g* A obtained within the two schemes differ by more than a factor of 2 
but the results for the observables calculated consistently within the two different 
approaches are typically within 1 % of each other. To be noticed that the presented 
SNPA two-nucleon weak current is not the only model available in the literature. In 
fact, in Ref. H3J two-body MEC's are derived from the hard pion chiral Lagrangians 
of the NAirpcoai system, and are not constrained to reproduce any experimental 
observable, like GT EXP . This is typically responsible for large mo del- dependence in 
the results, as some of the coupling constants and cutoff parameters entering the 
axial current are poorly known. 

The two-body axial current operator in %EFT consists of two contributions: a 
one-pion exchange term and a two-nucleon contact-term. The explicit expressions 
for these terms can be found in Ref. [2TJ While the coupling constants which appear 
in the one-pion exchange term are fixed by 7rA data, the LEC which determines 
the strength of the contact-term (da) has been fixed by reproducing GT EXP . The 
values of cIr for A=500-^800 MeV are given in Table [4]L39J experimental error 
on GT EXP is primarily responsible for the uncertainty in <1r. 

Our results of Ref. [39] for reaction ([7]) are compared in Table [2] with those of 
two previous calculations, performed in SNPA and xEFT*E3. The first one uses 
the Nijmegen I and Nijmegen 93 ^21 Hamiltonian models to obtain the nuclear wave 
functions, and MEC's derived from the Lagrangians of the NAnpuai system. The 
second calculation uses the AVIS^ potential to derive the wave functions, and the 
same X^FT weak current model presented above, constrained to reproduce GT EXP 
in tritium /3-decay. However, the l7rC, 27T and contact-term contributions to the 
weak vector current are not included. Furthermore, only the S-wave contribution in 
the nn final scattering state (the 1 Sq state) is retained, and higher partial- wave con- 
tributions are estimated based on Ref. [63l By inspection of Table [2] we can conclude 
that: (i) our calculated T D values are in good agreement with the results of Ando 
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et al. S3, and the small existing differences has been traced back to the inclusion 
in the weak vector current of the IttC, 2tt and contact-term contributions 1221 (jj) 
The calculated T D value of Ricci et al. differs from the other results by 7-f-lO 
%. In order to investigate the origin of the discrepancies between our results ^21 
and those of Ando et al. S3 on one side, and the results of Ricci et al. ^31 on the 
other, we have repeated ^21 the calculation of T D (and To), including the so-called 
"potential currents" , i.e. those operators arising when PCAC is implemented at the 
two-body level. It was argued in fact by Ricci et al. ^31 that "omitting the potential 
current causes an enhancement of the doublet transition rate Ai/ 2 [i.e. T D ] by ~ 1 
%" . These currents were first constructed in Ref. EH and we have recently reviewed 
them in Ref. E21 where their explicit expression can be found. The calculation has 
been performed within xEFT*, using the AV18 (AV18/UIX for A = 3) Hamil- 
tonian model. Again the LEC which determines the strength of the axial current 
contact-term has been fixed by reproducing GT EXP . The results for T D (and To) 
are 393.2 ± 0.8 s^ 1 (1488 ± 9 s^ 1 ), with r £, ( 1 5 )=250.1 ± 0.8 s" 1 , in perfect 
agreement with our previous results of Ref. 39. From this we can conclude that the 
potential currents proposed by Ricci et al. give negligible contributions to the 
rate T D (and To), and the discrepancy between the theoretical calculations is still a 
puzzling problem. From a historical point of view, it should be noticed that such a 
discrepa ncy betw een different theoretical resul ts fo r T D already existed in the early 
nineties ISSMIS] In fact; Adam and T^nk 66 found t d = 416 ± 7 s -i ) wm l e 

Tatara et al.^ and Doi et al. found T D = 399 s" 1 and 402 s _1 , respectively. 

Finally, we should also mention that a calculation of T D has been performed 
within pionless EFT by Chen et al. S3. The objective of this work, however, is not 
to predict T D , but rather to find the relation between the two-nucleon axial current 
matrix element entering the muon capture rate on deuteron and the pp weak capture. 
Within this approach, therefore, a precise experimental determination of T D will 
put a stringent constraint on this matrix element, and consequently on the pp weak 
capture rate. 

Using A — 3 nuclear wave functions derived from the AV18/UIX or N3LO/N2LO 
Hamiltonian models, and the same SNPA or xEFT* weak charge and current op- 
erators pres ented above, we have studied also the total capture rate Tq for reac- 
tion (|8]) 1221 The results are shown in Ta ble |3l and are compared with other theo- 
retical works of the last ten years I 40 | 45 | Q ur ca l cu lation of Ref. |40] represents the 
first attempt to study muon capture on 3 He in a way that is consiste nt wit h the 
approach adopted for the weak proton capture reactions pp and hep 1 1 9 | 20 |_ 
nuclear wave functions were obtained, within the HH method, from the AV18/UIX 
Hamiltonian model, and the nuclear weak current was derived within the SNPA, 
as presented above. The theoretical uncertainty reported in Table [3] for To results 
from the adopted fitting procedure and experimental err or on GT EXP . Note that 
a calculation based on the older Argonne U14 (AV14) ^3 two-nucleon and Tucson- 
Melbourne (TM) ESI Hamiltonian model yielded a r of 1486 ± 8 s -1 , suggesting a 
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weak model-dependence. In fact, we have demonstrated SOI that Tq roughly scales 
as the triton binding energy. Therefore, any meaningful comparison between results 
obtained using different Hamiltonian models requires the inclusion of three-nucleon 
forces. This is the reason why we have not considered in Table [3] the results of Ho 
et al. obtained, within the SNPA, without the inclusion of MEC's and, most 
important, three-nucleon interaction. In Ref. |40]we provide also the only available 
recent theoretical prediction for the recoil asymmetry A v , found to be 0.5350 ± 
0.0014 with the AV18/UIX, in agreement with the experimental result of Ref. |36j 
0.63 ± 0.09 (stat.)+°;^ (syst.). The results for A v are very little model-dependent, 
but very sensitive to Gps(q 2 )- A v would vary by roughly 20 %, if Gps{q 2 ) would 
be 50 % larger than the PCAC value (see Fig. 1 of Ref. |40| . The corresponding 
variation for To would be of the order of 5 %. 

The first study of reaction ([8]) within xEFT* approach has been performed by 
Gazit S3. The nuclear wave functions have been obtained with the Effective Inter- 
action HH method and the xEFT weak current presented above. However, as in 
Ref. 1331 no l7rC, 2ir and contact-term contributions to the weak vector current are 
retained. The theoretical uncertainty reported in Table [3] has two main sources: the 
experimental uncert aint y on the triton half-life, and the calculation of electroweak 
radiative corrections ^H. Few comments are here in order: (i) electroweak radiative 
corrections were not included in our studies of Refs. [3J2 and HOI • Were to be in- 
cluded, the central value for Tq would become 1493 s _1 for both calculations, in 
nice agreement with the result of GazitSS. (ii) The comparison between our study 
of Ref. |40] and that of Gazit ^ suggests that the SNPA and xEFT* results nicely 
agree, when the MEC's are constrained to reproduce GT EXP . We have verified this 
observation^^ for both reaction (J7]) and ([8]). However, we have shown that lirC, 2n, 
and contact terms in the mesonic %EFT vector current are important in orde r to 
achieve such an agreement. If they were to be neglected, Tq would be 1453 s _1 ^21. 

Finally, we recall the studies of the early nineties by Congleton and Fearing 
and Congleton and TruhlikL^U In the latter work, the nuclear wave functions were 
obtained from the AV14/TM Hamiltonian model and the nuclear weak current 
retained contributions similar to those of Ref. |44] The value obtained for the total 
capture rate To was 1502 ± 32 s _1 , the uncertainty due to poor knowledge of 
coupling constants and cutoff parameters. 

Only very recently, the first steps to study muon capture reactions in a con- 
sistent xEFT framework have been done I 48 | 49 |_ Although the results are not yet 
published, the main ingredients of a xEFT calculation are outlined. In particular, 
in Ref. [48j we have used the N3LO and N3LO/N2LO interaction models and the 
xEFT weak current operator presented above. Furthermore, the LEC cIr determin- 
ing the strength of the axial current contact-term, and the LEC cd, entering the 
contact-term three-nucleon interaction at N2LO, have been related, as suggested in 
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Refs.[72]and[73l as 



d R = -^^cd + \m n {c z + 2c 4 ) + ~ , (35) 
A x g A 3 6 



where C3 and C4 are the LEC's of the irN Lagrangian, already part of the chiral two- 
nucleon potential at NLO, and A x = 700 MeV is the the chiral-symmetry-breaking 
scale. Then, the calculation is implemented in the following steps: (i) all the LEC's 
present in the interaction and in the current are set consistently, and the same cutoff 
regulator, i.e. exp[— (q/A) 4 ], is used both in the current and in the N2LO three- 
nucleon interaction, (ii) The 3 H and 3 He ground state wave functions are calculated, 
within the HH method, u sing th e N3LO/N2LO Hamiltonian model, for two values 
of A = 500 and 600 MeV I ' \ and the set of values {cd,ce} are determined, for 
which the A = 3 experimental binding energies (BE's) are reproduced. A wide range 
of cd values has been spanned, and in correspondence to each cd in this range, Ce 
has been fixed in order to reproduce either BE( 3 H) or BE( 3 He). (iii) For each set 
of {cd,ce}, the triton and 3 He wave functions are calculated and, using the xEFT 
axial weak current discussed above, the Gamow- Teller matrix element of tritium 
/3-decay (GT TH ) is determined. This allows to determine the range of cd values 

rpTJ T7 1 Y F* 

for which GT = GT within the experimental error. A corresponding range 
for ce is given from the previous step, (iv) For the minimum and maximum values 
of {cd;ce} in the selected range, the LEC's g^s and g^y entering the two-nucleon 
contact terms of the electromagnetic current, and therefore the weak vector current, 
are determined by reproducing the A = 3 magnetic moments. At this point, the 
potential and current models are fully constrained, and the results for T D and Tq 
are xEFT predictions. They are found to be T D = 399 ± 3 s" 1 (r I? ( 1 S' ) = 255 ± 1 
s _1 ) and Tq = 1494 ± 21 s _1 , including electroweak radiative corrections These 
results are in good agreement with the ones of the other calculations mentioned 
above, except for that of Ricci et al. for which the discrepancy remains of the 
order of 44-9 %. On the other hand, in a similar calculation, Adam et al. ^2 have 
found T D in the range 401.2 s _1 4436.6 s _1 , depending on the xEFT two-nucleon 
potential used. 

We conclude remarking that a comparison between the calculated and measured 
rates for muon capture on 3 He makes it possible to put a constraint on the induced 
pseudoscalar form factor Gps{q 2 ) at Qq = — 0.954 mi, relevant for this reaction. 
A similar comparison could be done for the muon capture on deuteron. However, 
being the available experimental data so uncertain, such a comparison would be 
less significant. Within the xEFT approach, we have varied Gps (<?o) to match the 
theoretical upper (lower) value with the experimental lower (upper) value for Tq S3. 
This has allowed to obtain for Gps(qo) = 8.2 ± 0.7, in very good agreement with 
the xPT prediction of Eq. ©, which gives G™(<?o) = 7.99 ± 0.20 ^ 
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4. Conclusions 

Muon capture reactions on light nuclei, in particular deuteron and 3 He, have demon- 
strated to be an interesting, fruitful and controversial field of research, both exper- 
imentally and theoretically. The work on this subject has been extensive, and the 
last few years have seen even a growth of interest and research. At this point, the 
experimental situation can be summarized as follows: (i) the total rate for muon 
capture on 3 He, reaction ©, is very well determined, with an accuracy of 0.3 %, 
hard to be reached by any present theoretical calculation, (ii) The angular correla- 
tion parameters, or so-called recoil asymmetries, are poorly known. Only the vector 
asymmetry A v has been measured, but the experimental error is still very large, of 
the order of ~ 30 %. (hi) The rate for muon capture on deuteron, reaction ([7]), from 
the doublet hyperhne state is also poorly known, with experimental values which 
agree among each other, but have uncertainties of 6-j-10 %. Howev er, the ongoing 
experiment performed by the MuSun Collaboration at PSI l ^ l ^ l w jjj clarify the 
situation and determine the rate with a precision of 1.5 %. 

The theoretical situation is evolving very fast, and a large effort has been put 
in the past few years to reduce as much as possible the the oretica l uncertainty on 
the calculated observables. In particular, it has been shown I 39|40 | ^ na |- f- ne moc j e i_ 
dependence, relative to the adopted models for the nuclear interaction and weak 
currents, can be strongly reduced by fitting the unknown parameters of the nu- 
clear currents to some significant observables, as the tritium half-life and the A = 3 
magnetic moments. A crucial role in these calculations is played by the numeri- 
cal techniques used to calculate the few-body wave functions, with the considered 
accurate (and highly complex) Hamiltonian models. Without such a fundamental 
ingredient, all the calculations mentioned above would be affected by a much larger 
uncertainty. However, significant discrepancies remain between the available the- 
oretical calculations of the rate for muon capture on deuteron, reaction 0, and 
their origin is still to be understood. Finally, the first st eps toward a xEFT, and 
ultimately QCD-based, prediction have been made HS|49| 

Considering the two main motivations to study muon capture on light nuclei, i.e. 
(i) to provide significant tests for the theoretical frameworks used in the study of 
reactions of astrophysical interest not accessible experimentally, and (ii) to extract 
the value for the induced pseudoscalar form factor Gps(q 2 ) and ultimately validate 
the xPT predictions, further investigations are highly recommended. In particular, 
there are very few studies on the muon capture rea ctions on 3 He in t he two - and 
three-body breakup channels, both theoretically and experimentally E ^ | 77 | p ur _ 
thermore, an accurate measurement for the angular correlation parameters of muon 
capture on 3 He could put even a more stringent constraint on Gps(q 2 )- 
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